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If the dark matter in the universe is a self-gravitating Bose-Einstein condensate (BEC) with quartic 
self-interaction described by the Gross-Pitaevskii-Poisson system, the adhesion model, the Burgers 
equation and the cosmological Kardar-Parisi-Zhang (KPZ) equation that have been introduced 
heuristically to solve the problems inherent to cold dark matter (CDM) models find a natural 
justification and an interesting generalization. 
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I. INTRODUCTION 



The large-scale structure of the universe can be rela- 
tively well-explained by cold dark matter (CDM) models. 
In these models, the dark matter is represented by a colli- 
sionless gas without pressure (dust) interacting via New- 
tonian gravity only. It can be described by the Vlasov- 
Poisson system or, under additional approximations, by 
the pressureless Eulcr-Poisson system in an expanding 
background However, this model faces several prob- 
lems on the scale of galactic or sub-galactic structures. 
For example, CDM simulations lead to dark matter ha- 
los with density cusps [H while observations of rotation 
curves favor flat density profiles. In addition, it predicts 
an abundance of satellite galaxies around each galactic 
halo that is far beyond what we see around the Milky 
Way. At the cosmological level, nonlinear gravitational 
clustering is often studied by resorting to the Zeldovich 
approximation [Sj. This approximation leads to the in- 
viscid Burgers equation which describes the free motion 
of fluid particles However, this equation develops 
shocks and formal singularities known as caustics (hav- 
ing the form of pancakes) , so that it becomes invalid after 
particle-crossing. If we accept multi-stream solutions, the 
particles just cross each other after the caustic and pass 
through the pancakes instead of clustering into smaller 
objects like groups of galaxies. Therefore, a small-scale 
regularization must be introduced to overcome this prob- 
lem and model the effects of "punctuated" gravitational 
attraction and pressure gradients that are not captured 
by the Zeldovich approximation. Gurbatov et al. Q in- 
troduced the so-called "adhesion model" in which parti- 
cles move according to the Zeldovich approximation until 
they fall into pancakes when their trajectories intersect, 
then "stick" to each other. This sticking can be mod- 
eled by introducing a viscosity in the Burgers equation 
in order to represent the effect of strong gravitational 
forces or pressure effects in the vicinity of a caustic. Of 
course, the viscosity must be small in order to provide 
a smoothing effect at small-scales only (where particle- 
crossing occurs) but the limit — >■ is different from 
taking v = Q. This model gives very good results in the 
nonlinear regime and can reproduce the skeleton of the 
"cosmic web" of large-scale structures (sheets, filaments, 
nodes) in iV-body numerical simulations (see, e.g., Ref. 
0). A completely different approach was developed by 



Widrow & Kaiser who proposed to describe a classi- 
cal coUisionless self-gravitating gas by the Schrodinger- 
Poisson system. In this approach, the constant Ti is not 
the Planck constant, but rather an adjustable parame- 
ter that controls the spatial resolution \deB through a de 
Broglie relation A^eS = h/mv. It is argued that when 
?i 0, the Vlasov-Poisson system is recovered and that 
a finite value of fi provides a small-scale regularization 
of the dynamics. In that case, the Schrodinger-Poisson 
system has nothing to do with quantum mechanics since 
it aims at describing the evolution of classical collision- 
less matter under the influence of gravity (in static or ex- 
panding universes) . This model was further developed by 
Short & Coles Q who introduced a free-particle approxi- 
mation. They showed that the dynamics of the particles 
before crossing is relatively close to the Zeldovich approx- 
imation but when crossing occurs the quantum pressure 
prevents the formation of singularities. Therefore, the 
quantum pressure replaces the role of the viscosity in the 
adhesion model. Although these models have given inter- 
esting results and can be very useful in practice, one can 
argue, however, that their justification remains relatively 
ad hoc. 

Recently, Bohmer & Harko @ have proposed that dark 
matter could be a self-gravitating Bose-Einstein conden- 
sate (BEC) with quartic self-interaction described by the 
Gross-Pitaevskii-Poisson (GPP) system. This idea takes 
its ori gin in the concept of boson stars introduced by 
Kaup |To| and Ruffini & Bonazzola [ll|. It is well- 
known that the Gross-Pitaevskii equation (or the non- 
linear Schrodinger equation) can be reduced to hydrody- 
namical equations by means of the Madelung jl2| trans- 
formation. This yields the Euler-Poisson system with a 
"classical" isotropic pressure due to self-interaction and 
a quantum anisotropic pressure arising from the Heisen- 
berg uncertainty principle. For a quartic self-interaction, 
the barotropic equation of state is that of a polytrope 
of index n — 1. At large-scales, the pressure terms are 
negligible and one recovers the CDM model which has 
proven very successful. However, at small-scales, the 
classical and quantum pressures can play a crucial role 
and regularize the dynamics by preventing the formation 
of singularities (caustics). BEC dark matter has very 
interesting properties: (i) the pressure can stabilize the 
system against gravitational collapse and lead to dark 
matter halos with flat density profiles. For very light 
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bosons without self-interaction, the stabilization is due 
to the quantum pressure. For more massive bosons with 
self-interaction, the stabilization is due to repulsive scat- 
tering. Therefore, a BEC dark matter can solve the cusp 
problem and the missing satellite problem (ii) 
BEC dark matter halos can reproduce the rotation curves 
of several low surface brightness galaxies [l5| . (iii) At a 
cosmological level, a BEC dark matter can accelerate the 
formation of structures with respect to ordinary CDM 
models [l^, In this brief report, we make the ad- 
ditional remark that if dark matter is a self-gravitating 
BEC with quartic self-interaction, the adhesion approxi- 
mation, the Burgers equation and the cosmological KPZ 
equation that have been introduced phenomenologically 
to solve the problems inherent to CDM models find a 
natural justification and an interesting generalization. 



II. STOCHASTIC 
GROSS-PITAEVSKII-POISSON SYSTEM 

We assume that dark matter is a self-gravitating BEC 
with quartic self-interaction described by the stochastic 
Gross-Pitaevskii-Poisson (GPP) system 
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where ip{r,t) is the wave function, p{r,t) = Nm\ip\'^ 
the density, $(r,t) the gravitational potential, g = 
4:TTash^/m^ the pseudo-potential accounting for short- 
range interactions (a^ is the s-scattering length) [l^ and 
r]{r,t) a stochastic potential (noise). We write the wave 
function in the form ^(r,i) = A{r,t)e^^^'^'*''^^ where A 
and S are real, and make the Madelung [l^ transforma- 
tion p = A^to|-0P = NmA^ and u = WS/m, where p{r, t) 
is the density field and u(r, t) the velocity field. We note 
that the flow is irrotational since V x u = 0. With this 
transformation, the stochastic GP equation ([!]) is equiv- 
alent to the stochastic barotropic Euler equations with 
an additional term Q = — (?i^/2m)A.yp/yp called the 
quantum potential (or quantum pressure). Indeed, one 
obtains the set of equations 
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The first term on the r.h.s. of the Euler equation Q can 
be interpreted as a classical isotropic pressure — (l/p)Vp 
described by a barotropic equation of state p = p{p)- For 
a quartic self-interaction, it is given by 
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This is a polytropic equation of state of the form p = Kp'^ 
with polytropic index n = 1 (i.e. 7 = l-|-l/n = 2) and 
polytropic constant K = g /2 = 2-Kasfi' jw? . Other types 
of barotropic equations of state can be obtained depend- 
ing on the form of the self-interaction. A detailed study 
of the time-independent solutions of the GPP system (or 
quantum barotropic Euler-Poisson system) , connecting 
the non-interacting limit (os ~ 0) [ll|] to the Thomas- 
Fcrmi TTF) limit (Q ~ 0) [| has been given recently in 
Ref. [3. 



III. BEC EQUATIONS IN AN EXPANDING 
UNIVERSE 



For the sake of simplicity, we consider an expandin 
Einstein-de Sitter (EdS) universe (A = and k = 0) [l| 
described by the equations 



Pha 



1, 



-TrGpfefl, b? = -TrGpba'^, (6) 



yielding a ex t^/s^ H ^ a/a = 2/3t and pb = l/GirGt'^, 
where pb{t) is the background density, a{t) the scale fac- 
tor and H = d/a the Hubble constant. Our approach can 
be easily extended to more general models of universe. 

We shall first rewrite the hydrodynamic equations in 
the comoving frame. Making the change of variables 
r = a(t)x, u — Hr + v and (j) = $ — where v is the pe- 
culiar velocity and 3>h = (2/3)7rG/5&(t)r^ the background 
gravitational potential, and introducing the density con- 
trast 5(x, t) = [p(x, t)~ pb{t)]/ pb{t), we obtain the system 
of equations [131 ■ 

— + -V • [(1 + d)v] =0, A(^ = 4TTGpba^6, (7) 
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For simplicity, we shall always denote the noise by 77 al- 
though a new notation should be introduced after each 
transformation. At large scales, pressure and noise ef- 
fects are negligible and the CDM model {h ^ p = 1] = 0) 
is recovered. However, pressure effects become impor- 
tant when nonlinear structures form, and the BEC dark 
matter model 0-® should be used instead. 

Measuring the evolution in terms of a rather than in 
terms of t and introducing the new velocity w = v/ad 
and the new gravitational potential = (j> / AttG pba^ , we 
can recast the foregoing equations in the form 
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For short times, the perturbations are small S <^ 1, 
(/) ^ 1, |w| ^ 1, and the hydrodynamic equations can be 
linearized. Wo obtain 
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Taking the "time" derivative of Eq. pH-a) and the di- 
vergence of Eq. (fT2l). these equations can be combined 
into a single equation for the density contrast 
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This equation has been studied in Ref. [17[ in the no- 
noise limit. At large scales, we can ignore pressure and 
noise effects, and wc recover the equation for the density 
contrast of a cold gas 
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The growing solution to this equation is 5+ (x, a) = 
aD{x) Q. Then, Eq. ((TTJb) implies that tp{^,a) = 'ip{x) 
is constant. Therefore, (5+(x, a) = aAip{x.). On the other 
hand, in the cold gas approximation, Eq. (jl2p reduces to 
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After a transient regime, the velocity field tends toward 
the solution w(x, a) ~ — VV'(x). 



IV. ZELDOVICH APPROXIMATION AND A 
GENERALIZED EQUATION OF STRUCTURE 
FORMATION 

The Zcldovich approximation Q amounts to extend- 
ing this relation to the (weakly) nonlinear regime, i.e. 
w(x, a) ~ — V'0(x, a). The fact that w is a potential 
flow in the BEG model gives further support to this ap- 
proximation. In that case, Eq. (fTO|) reduces to the form 



9w 



3r 



(w • V)w = — 



Sash. 
2Gm3a4 



V(5 



^fAVTTs\ ^ 

V I . I - V?7, 



WTTGm'^Pba^ V VTTS 



(16) 



where the explicit dependence on the gravitational po- 
tential has disappeared. Since S = aAijj = — aV • w and 



V5 = — aAw, we can rewrite the foregoing equation as 
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This can be viewed as a generalized noisy Burgers equa- 
tion with an additional quantum pressure term. Since w 
is a potential flow, we obtain 
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This can be viewed as a cosmological KPZ equation [l9j . 
Without the forcing term, this is just the Bernouilli equa- 
tion, or the Hamilton- Jacobi equation, with an additional 
quantum potential. We note that the potential "0 is re- 
lated to the phase S of the wave function of the BEG by 
-0 = —{1/ a){S/ma^H — /2). A noisy Burgers equation 
and a KPZ equation have been introduced and studied 
previously in cosmology by Buchert et al. [201 . Jones 
[2l| . Goles [12 and Matarrese & Mohayaee [23l|. 
Let us consider particular cases of Eq. ([T7|). 

(i) For tts = h = r] = 0, we get the inviscid Burgers 
equation which is equivalent to the Zeldovich approxima- 
tion. 

(ii) If we neglect the quantum potential, we get the 
noisy Burgers equation 

dw 

— + (w • V)vir = :/(a)Aw - V?7, (19) 
da 

with a time-dependent viscosity given by J^(a) = 
iagh^ /2G'w?a^ (this equation is time-reversible via a ^• 
—a and w — > — w). Therefore, the BEG dark matter hy- 
pothesis leads to a natural justification of the (stochas- 
tic) adhesion model. Another justification was previ- 
ously given by Buchert et al. [20| in terms of a coarse- 
graining process inherent to a hydrodynamic descrip- 
tion. However, to arrive at the Burgers equation, several 
simplifying hypothesis had to be introduced: (i) dark 
matter is described by hydrodynamic equations with an 
isotropic pressure; (ii) the equation of state is a poly- 
trope p = Kp"' with 7 = 2; (iii) the velocity field and 
the stochastic force derive from a potential. Interest- 
ingly, these properties directly result from the BEG dark 
matter hypothesis (i.e. from the GP equation) without 
further assumption. Note also that the change of variable 
tp{x,a) = 2i'{a)\nW{x,a) transforms Eq. ([T8[) into 



— — = iy(a)AW- 
da 



^77(x,a) + -lnW^ 
2iy a 
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(iii) In the absence of self-interaction and noise, we 
obtain an equation of the form 

t2 



dw 

— h (w • V)w = 

da 



3r 
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This can be viewed as a Burgers equation where the 
small-scale regularization is provided by the quantum 
pressure. Our approach has some similarities with the 
effective wave mechanics approach of Short & Cole Q. 
However, in our approach, the coefficient in front of the 
quantum potential in Eq. ([2T|) depends on "time" a, 
while in the approach of Short & Cole [1], this coeffi- 
cient is constant (see their Eq. (24)). This is because 
they introduce an effective Schrodinger equation [their 
Eq. (23)] directly in the comoving frame while we start 
from the (nonlinear) Schrodinger equation ([1]) in the in- 
ertial frame. When we write the (nonlinear) Schrodinger 
equation ([!]) in the comoving frame, we obtain an equa- 
tion different from their Eq. (23). Ribeiro & Peixoto 
de Faria [IJ] have also developed an effective wave me- 
chanics approach in which they relate the gradient of 
the quantum potential to the Laplacian of a kinemati- 
cal velocity plus a noise term. Their approach, which is 
essentially phenomenological, leads to a result different 
fromEq. 

V. CONCLUSION 

In this brief report, we have shown that the assump- 
tion that dark matter is a self-gravitating Bose-Einstein 
condensate with quartic self-interaction described by the 
(stochastic) GPP system leads to a natural justification 
of the adhesion model, the Burgers equation and the cos- 
mological KPZ equation without more hypothesis than 
the Zcldovich approximation that is common to most 
works on the subject. Therefore, not only the BEC model 
is consistent with previous works, but it generalizes them 
and extends their scope. In addition, it gives a new jus- 
tification of these phenomenological models. This result 



adds to the other nice properties of BEC dark matter (flat 
density profiles, flat rotation curves, acceleration of the 
growth of perturbations...) recalled in the Introduction. 

The BEC model not only follows the general evolu- 
tion of inhomogcneities but it also describes the internal 
structure of density enhancements. In this sense, it im- 
proves upon the standard adhesion model. Indeed, at 
the level of dark matter halos, the equations reduce to 
the condition of hydrostatic equilibrium and lead to viri- 
alized structures similar to n = 1 polytropes (other con- 
figurations could be obtained depending on the form of 
self- interaction and on the equation of state). Therefore, 
BEC dark matter provides a model that describes both 
the large-scale structures of the universe (through the 
generalized Burgers equation ([TT])) and the structure of 
dark matter halos (through the explicit expression of the 
pressure ^ arising from short-range interactions). 

If we justify the GPP system in terms of quantum me- 
chanics, we must assume that the mass m of the bosons 
is extraordinarily small (to < 10~^^eV/c^!) for quantum 
mechanics to be relevant on galactic scales. This has been 
called "fuzzy dark matter" in Ref. [31 . Alternatively, 
one can produce similar quantum effects with larger bo- 
son masses if the particles have a self- interaction [25j . 
Since the nature of dark matter is not known, we can- 
not reject a priori that quantum mechanics plays some 
role at galactic or cosmological scales. In any case, the 
GPP system (HJ-lH]) can always be viewed as an effec- 
tive wave mechanics approach with an adjustable Planck 
constant which generalizes the initial model of Widrow & 
Kaiser j7j based on the Schrodinger-Poisson (SP) system. 
Finally, other approaches to the problem, based on the 
Vlasov-Poisson (VP) system, could be contemplated, see 
e.g. Ref. [13. 
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